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Abstract. - The spin-spin correlation length and the static structure factor for bilayer antifer- 
romagnets, such as YBa2Cu306, are calculated using field theoretical and numerical methods. 
It is shown that these quantities can be directly measured in neutron scattering experiments 
using energy integrated two-axis scan despite the strong intensity modulation perpendicular to 
the layers. Our calculations show that the correlation length of the bilayer antiferromagnet 
diverges considerably more rapidly, as the temperature tends to zero, than the correlation 
length of the corresponding single layer antiferromagnet typified by La2Cu04. This rapid 
divergence may have important consequences with respect to magnetic fluctuations of the doped 
superconductors. 



A powerful method to measure the spin-spin correlation length of a layered magnet is the 
neutron scattering method known as the energy integrated two-axis scan (TAS). In recent 
years, TAS has been successfully applied to L^CuO^ji]], which is the parent compound of 
one of the high temperature superconductors. This experimental technique is not readily 
extendable, however, to a wide class of high temperature superconductors with close magnetic 
bilayers or triple layers within the unit cell; a particularly important example is YBa2Cu30g, 
which has a close pair of magnetic planes within the unit cell, but it is otherwise a square 
lattice spin S — 1/2 Heisenberg antiferromagnet. The reason for the difficulty is an intensity 
modulation^ with the momentum transfer perpendicular to the planes. Thus, there have 
been no direct measurements of the correlation length despite considerable discussion of the 
importance of antiferromagnetic fluctuations in these materials. 

In the present paper, we use both the field theoretical approachQ and the numerical 
quantum Monte Carlo (QMC) loop algorithm^ to obtain the low temperature properties of 
bilayer antiferromagnets. We also show that in the experimentally relevant regime TAS can be 
extended to such antiferromagnets. Thus, it is hoped that the antiferromagnetic fluctuations 
can be explored more thoroughly in future measurements. This is likely to be important 
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in understanding the magnetic properties of the superconductors obtained by doping these 
antfcrromagnetic parent compounds. 

The Heisenberg model for a spin-S* bilayer antifcrromagnet is 



lip) 



(i) 



j(2) 



(1) 



The sum in the first term is over the nearest neighbor pairs on a square lattice in each plane, 
where the plane index p takes two values 1 and 2. The second term represents the coupling 
between the planes. The exchange constants J\\ and J± are both positive. 

The low energy, long wavelength properties of the two dimensional Heisenberg model is well- 
described by the quantum 0(3) nonlinear ct- model ||. Here, we shall consider its generalization 
to coupled bilayers. The Euclidean action for this system can be written down on general 
symmetry grounds^, but it can also be derived from a (1/5) expansion ||. The action is 



s = r d XQ [d 2 x v (j-\d^\ 2 -h *w 



(2) 



where is the component of the staggered field fl^ in the direction of the staggered 
magnetic field ho- Here, all the coupling constants and dimensional variables have been scaled 
to their dimensionless forms. We shall work directly at d = 2; the label \i = 1, 2 denotes the 
spatial directions, and fi = denotes the imaginary time direction of extent u. 

The relations between the er-model parameters and the Heisenberg model parameters are 
known in the large-S" limit ^ q], but we shall not need them in the present papper. The 
momentum cutoff of the a-model, A, is chosen to be V2n/a, where a is the lattice constant of 
the Heisenberg model, to conserve the number of degrees of freedom. The bare dimensional 
bilayer gap is given by Aj_ = TicoK.yj2goj ' cf^, where cq is the bare spin wave velocity at the 
scale A -1 . 

It can be seen from explicit calculations JtJ that the angular momentum coupling between 
the layers is irrelevant for weakly coupled bilayer systems. Moreover, this coupling reduces to a 
higher gradient coupling when the number of layers tends to infinity. Therefore, we shall omit 
this term even though its presence breaks the "Lorentz invariace" and slightly renormalizes 
the spin-wave velocity. 

For notational simplicity, it is useful to define ho = hogo and 7^ = (2go)/gj_. The one-loop 
momentum-shell calculations similar to those of Ref. p| yields 
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2/j, (6) 
where e l is the length rescaling factor, F 1 (g/2t,h) = coth^V 1 + h/2t)/\/l + h, and F 2 (g/2t, h, j±) 



2j±-^F 2 (g/2t,h,j x ), 



coth(py 1 + h + J±/2t)/\/ 1 + h + j±. The variable t is the dimensionless temperature vari- 
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able. The dimcnsionlcss thickness of the slab in the imaginary time direction u = g/t satisfies 
a simple scaling relation given by (g/t) = (go/to)e~ l . 

The zero temperature flows of the renormalization group equations are shown in Fig. 1. 
There are two phases, separated by a separatrix between the two unstable fixed points (g = 
4:TT,g± = oo) and (g — 8ir,g± = 0); the former is the fixed point of the single layer casej|. 
There are two stable fixed points. The ordered-phase fixed point is located at (0,0), where 
both the in- and the inter-plane couplings are infinitely strong. The disordered-phase fixed 
point is located at (oo,oo), where the system becomes totally disordered. Although quantum 
nonlinear cr-model is a very accurate description of the low energy physics in or near the 
ordered phase, far into the disordered phase, such a continuum theory can not be expected to 
be valid. 

From experiments on YBa2Cu306+ x , it is known that the ground state is an ordered Neel 
state and that Jj_ ~ 0.1J||||]. Thus, the parameters are such that 7^ ~ (Jl/Jii) <C 1, and 
g is well below the critical value required for the phase transition to the quantum disordered 
phase at T = Q. Therefore, the system is in the renormalized classical regime 

To proceed further, we need an analytical solution of the renormalization group equations. 
We have obtained a good appoximation to the solution based on the following observations. 
In the renormalized classical regime, the bilayer gap VxT is initially much smaller than unity, 



but increases as 



7_l oc e 



a(l)l/2 



, where a(l) < 2, but tends to 2 for large We can therefore 
consider two regions, 71 <C 1 in region (I), and ~ 277. in the region (II). We solve 
the renormalization group equations separately in regions (I) and (II), and then join the two 
solutions together to get the final answer. The result is 



t t 4ir 



sinh 2 (-^ 
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sinh(^e-<)sinh (f^e^ + 7° (1 - M)) 



(7) 



The isolated single layer result is trivially recovered when the bilayer gap is much smaller than 
the inverse correlation length, that is, \f^\_ <C er l ■ In the limit of large I, Eq. (fjj) has the 
asymptotic form 2/t = 2/t e g — l/2w, where 
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From the solution of the renormalization group equations at T — 0, the quantity inside the 
first square parenthesis on the right hand side of Eq. (||) can be shown to be the one-loop 
renormalized zero temperature bilayer spin stiffness constant. Similarly, one can show that 

Aj_(0) = hcK^j 7^(1 — f^) is the one-loop renormalized bilayer gap at T = 0. We assume that 
these renormalizations hold to all orders. Thus, we can express t e s in terms of the physical 
bilayer parameters, denoted by the superscript 6, as 
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t e fi 
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The problem is now mapped onto an effective classical cr-model if we identify 2/t c g = 1/tf 
and 2/t = l/t cl . One can then show[p| that 



-A — e*eff 
8 4tt 
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with the exact prefactor determined in Ref. 



4 



EUROPHYSICS LETTERS 



The physical parameters can be calculated from the spin wave theory of the antiferro- 
magnetic bilayer Heisenberg model in Eq. (|l|). For S = 1/2, the spin wave velocity of the 
bilayer complex calculated to order (1/5*) is given by He = V2J\\Z c a(l + jj~) 1/2 ; similarly, 

A b ± (0) = 2y/JplZ A . For J ± = 0.08J|| §, Z A = 0.95 and Z c = 1.15. The bilayer spin stiffness 
constant can be calculated from the hydrodynamical relational p b s (0) = ^J\\Z%Z X , where the 
renormalization factor for the uniform susceptibility Z x = 0.53 for the same set of parameters. 

The static structure factor contains two pieces, one for the symmetric combination of the 
unit vectors from each layer and the other for the antisymmetric combination. The notation 
in the cr-model is exactly the opposite to the notation in the Heisenberg model, since the unit 
vector field fl^ in the a- model is the continuum limit of the direction vector of the staggered 
spin operator (— l) t+p S^ in the Heisenberg model. For the symmetric piece, we find that 
S s (k,t ) = 4£ 2 f(x)/(2ir), where x = fc£ and f(x) = [1 + |ln(l + a; 2 )](l + a; 2 )" 1 . Similarly, 
the antisymmetric piece is S a (fc,io) — t (2 — ^){k 2 + ^J 2 ) -1 . where £j_ is the length scale 
associated with the bilayer gap; approximately, we have £j_ = 7ic/Aj_(0). The symmetric piece 
is clearly dominant in the long- wavelength limit. 

In TAS, the wavevector of the incoming neutron is fixed, while the outgoing neutrons in a 
direction perpendicular to the layers are collected, regardless of their energies. The transferred 
wavevector is given by q = qy — q^. Its in-plane component is a constant, q|| = — q^ii, while 
its perpendicular component is a variable, q± = qf — qi±. For near the reciprocal lattice 
vectors, the form factor is approximately a constant, and the intensity is an integral over 
the 3D dynamic structure factor, which is related to the 2D dynamic structure factors by 
^)5 2D ( qj , W )+COS 2 (^)5 2D (q||, 

[lb. 

respect to the layers, symmetric in the a- model sense, and S 2 (qii,a;) to the corresponding 
symmetric spin combination, antisymmetric in the cr-model sense. 

In experiments one probes the region « G, where G is the nearest antiferromagnetic 
reciprocal lattice vector. Defining k = q|| — G, we can rewrite the intensity in terms of the 
cr-model structure factors. The intensity consists of two parts, I(k) = / s (k) + 2a. (k), where 



S* 3D (q, lj) = sin 2 (^)5 2D (qj|, to) + cos 2 (2^)S' 2D (q|| , to), where h is the distance between the 
two layers. The quantity S 2 (qii,o;) corresponds to the antisymmetric spin combination with 



J 8 (k) ~ / e%sm z 
Jo 

j* oo 

J a (k) ~ / dq f cos 
Jo 



2 

0/ - H-L)h 



^fr^tf "«?)), (11) 

Sr°(k,^(q} -<£)). (12) 



The qx-modulation is unimportant in the critical region. The reason is that S 2T)<7 (k, oj) 
is dominated by the critical fluctuations near uj = 0, where both sin 2 ^ Qf ^-O^ and 

, fc2 

2^ = j^q/ are essentially constants. Therefore, we can pull these factors out of the 
integal and obtain the intensity approximately proportional to the static structure factor, 
J 8 (k) ~ ^-sin 2 ( fe 7'- L) /i) f™ Et duS* D "Qt,w) ~ S s (k,t ), where we have reverted to the 
previous notation by dropping the superscripts. The quantity Ei is the incident neutron 
energy. For the antisymetric piece, we get an upper bound by neglecting the factor cos 2 (^— ), 
which is dqfSl D<T (k, ^(q 2 - q 2 )) ~ ^ ^ 5 a (k, to), where (qf) is some average of the 
wavevector. Because S s (0,t ) 3> 5 a (0,io) ; the intensity is dominated by the contribution 
from the symmetric piece. Therefore TAS for a bilayer should yield information about the 
symmetric piece, hence the correlation length. The contribution of the antisymmetric piece 
should result in a small broad background. 

The NMR relaxation rate for the in-plane Cu site will be given by[[To) ^ cx Ta£, because 
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Fig. 1. - Left: Zero temperature flow diagram. 



Fig. 2. - Right: Correlation lengths of the bilayer antiferromagnet compared to those of the isolated 
single layer. The solid circles are the Monte Carlo results. 



the symmetric structure factor will dominate when the correlation length £ is large. Thus, the 
formula for the correlation length can also be tested in NMR experiments, as for the single 
layer La 2 Cu0 4 @. 

For single layer antiferromagnets, QMC simulations have shown that the field theoretical 
expression for the correlation length is accurate for only very large correlation lengths, often 
larger than those accessible in experiments |l2| . To obtain results that are also reliable at 
intermediate and high temperatures we use QMC j| to calculate the correlation length defined 
from the second moment of the structure factor. As a function of the system size, the results 
converge within our statistical errors for systems of linear dimension L larger than 7£. We 
take this into account in estimating the infinite volume correlation length. 

The definition of the correlation length from the second moment of the static structure 
factor is not equivalent to that used in the field theoretical approach. The reason is the 
logarithmic correction 1 + ■kBf ln(l + x 2 ) to the Lorentzian form, where Bf = 1 in one loop 
order. However, from previous simulations, it is known that Br is strongly renormalized and 



that the actual correction is an order of magnitude smaller) 13 . The two definitions of the 
correlation length thus differ by only a few percent. 

In Table 1 we present our results for an inter layer coupling of J± = 0.08J|| and plot 
them in Fig. 2 together with the field theoretical results. Not unexpectedly there are 



Table I. - Second moment correlation length of the bilayer antiferromagnet at an inter layer coupling 
J_l = 0.08J|| calculated in quantum Monte Carlo simulations. 



T/J\\ 0.35 0.36 0.38 0.4 0.45 0.5 0.6 0.7 1 

£ 59.7(4) 43.3(3) 24.5(2) 15.6(1) 7.14(7) 4.37(4) 2.46(2) 1.73(1) 0.96(1) 



differences at intermediate temperatures, but the correlation length rapidly approaches Eq. 
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( [To| ) at lower temperatures. We find it surprising that at temperatures as high as 5Jj_(= 0.4J||) 
the correlation length of the bilayer antiferromagnet is already significantly higher than the 
correlation length of the single layer antiferromagnet [fl4| . 

In summary, the theoretical results obtained here should be important in understanding the 
magnetic fluctuations of bilayer high temperature superconductors, as they can be explored in 
neutron scattering experiments using the energy integrated two axis scan. The applicability 
of this technique is due to the fact that the spin-spin correlation length at low temperatures is 
exceedingly long in bilayer antiferromagnets[]l5|, and nearly critical fluctuations dominate for 
the parameter regime that is of experimental interest. The energy integration is thus correctly 
carried out. The present work shows that the bilayer complex in high temperatures such as 
YBCO are effectively strongly coupled despite the fact that J± ~ O.lJy. 



We thank O. Syljuasen for interesting discussions. This work was supported by the National 
Science Foundation, Grant. No. DMR-9531575. The QMC calculations were performed on 
the Hitachi SR2201 massively parallel computer of the computer center at the University of 
Tokyo. 
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